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. ABSTRACT 

O ■ 

, Aims. We have analysed ground-based multi-colour Stromgren photometry and single-filter photometry from the star tracker on the wire 
O |' satellite of the S Scuti star e Cephei. 

Methods. The ground-based data set consists of 16 nights of data collected over 164 days, while the satellite data are nearly continuous coverage 
of the star during 14 days. The spectral window and noise level of the satellite data are superior to the ground-based data and this data set is 
used to locate the frequencies. However, we can use the ground-based data to improve the accuracy of the frequencies due to the much longer 
time baseline. 

Results. We detect 26 oscillation frequencies in the wire data set, but only some of these can be seen clearly in the ground-based data. We have 
. ' usec * ' ne multi-colour ground-based photometry to determine amplitude and phase differences in the Stromgren b — y colour and the y filter in 
K> an attempt to identify the radial degree of the oscillation frequencies. We conclude that the accuracies of the amplitudes and phases are not 
» I ■ sufficient to constrain theoretical models of e Cep. We find no evidence for rotational splitting or the large separation among the frequencies 
1 detected in the wire data set. 

Conclusions. To be able to identify oscillation frequencies in 6 Scuti stars with the method we have applied, it is crucial to obtain more 
complete coverage from multi-site campaigns with a long time baseline and in multiple filters. This is important when planning photometric 
and spectroscopic ground-based support for future satellite missions like corot and kepler. 

Key words. Stars: oscillations; Stars: variable: <5 Set; Stars: individual: e Cephei (HD 211336; HR 8494) 



1. Introduction 

6 Scuti stars are main sequence population I A- and F-type 
stars. They are found in the classical Cepheid instability strip 
on the main sequence, have masses around twice solar, and 
temperatures around 7 500 K. They have only very shallow 
outer convection zones but their cores are fully convective. 
Most stars found in the instability strip are variable and most 
are multi-periodic, with periods around 1-2 hours and ampli- 
tudes typically at the mmag level, alth ough some have ampli- 
tudes above 0.1 mag (Rodri guez et all EoOO). 

To get a better understanding of 5 Scuti stars, several multi- 
site ground-based campaigns have carried out extensive moni- 
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toring of selec ted targets. F or exa mple, the Delta Scuti Network 
(Zima, 1997; Zim a et all 12002) has monitored FG Virginis 
duri ng severa l seas ons and 79 frequencies have been detected 
jBreeer et all 120051) . In that case, data from several ground- 
based observatories collected over 13 years were combined, 
which allowed the detection of frequencies down to amplitudes 
of just 0.2 mmag in Stromgren y. It is important to note that 
about 50 of th e reported fre quencies have amplitudes below 0.5 
mmag. IBreger et all J2005) stated that typical multi-site cam- 
paigns with a duration of 200 to 300 h detect only 5-10 fre- 
quencies. Since so few frequencies are detected, while t heory 
predic ts a much higher number of excited modes. IBreger et all 
(2005) concluded that either longer ground-based photometric 
campaigns or high-precision space-based campaigns are neces- 
sary. 
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Fig. 1. Each panel shows the light curve of e Cep from wire and osn in the uvby filters and b - y. The top panel shows the 
complete set of light curves and the bottom panel shows the details of a small part. Note that the osn data were collected about 
2.9 years prior to the wire data, i.e. the zero points fo are different (see main text). 



Buz asi et"afl d2005l) recently found Altair to be a low- 
amplitude 6 Scuti star based on high-precision photometry 
from the star tracker on the Wide-field InfraRed Explorer 
(wire) satellite. They found only seven frequencies but three 
of these had amplitudes around 0.1 mmag. The reason why so 
few frequencies were detected in Altair compared to FG Vir 
may be the much higher rotation rate of Altair, which has 
vsim = 210 ± 20 kms" 1 while FG Vir has vsin; aro und 
21 km s" 1 dMantegazza & Porettlll2002tlMittermaver & WeissL 
120031) . Suarez et al. (2005) have calculated models of Altair 
and used the cons t raints from the fr equencies detected by 
iBuzasi et all d2005l) . ISuarez et ail d2005l) found that high rota- 
tional velocity makes the interpretation of the frequencies diffi- 
cult, due to the limits of the applied second order perturbation 
theory and effects of near degeneracy. The 6 Scuti star treated 
in this paper has vsin; g 90 kms" 1 and m ay p resent a sim- 
pler ca se. The studies of iBreger et al.l d2005l) and lBuzasi et alJ 
(2005| both agree that high photometric precision and long 
temporal coverage is needed to fully explore the oscillation 
spectra of 5 Scuti stars. 

In the present study we have combined two quite different 
photometric data sets for the multi-mode 5 Scuti star e Cephei 
(HD 211336). One data set is from the wire satellite covering 
14 days with very high signal-to-noise (S/N). The other data set 
is single-site ground-based Stromgren uvby photometry from 
16 nights collected during 164 days, thus having gaps from a 



single day to several weeks. We find that by combining the very 
different properties of the data sets in terms of S/N and spec- 
tral window, we can measure the frequencies very accurately. 
We assess the uncertainties on the measured frequencies from 
extensive simulations to see if mode identification is possible 
using amplitude ratios and phase differences in the modes mea- 
sured with the Stromgren filters. 

1.1. Previous studies of e Cep 

e Cep was first identified as a 5 Scuti star by iBregerl dl966l) 
who found a single frequency at / = 23.8 ± 1.7 c/day based 
on only two n ights of data. One night of y band photometry 
bv lFesenl d 19731) showed an apparent change in amplitude, but 
today we know that this is likely due to beating of frequen- 
cies. Spectroscopy was first carried out during three nights by 
Gray| dl97lh . who confirmed the period known from photome- 
try. The mean radial velocity amplitude measured on the three 
nights was 15 + 3 kms -1 , but this disagrees with upper limits 
of ~ 1 km s _1 found bv lKennellv et alJ dl999l) and lBaade et alJ 
dl993l) . 

Line p r ofile variations (LPVs) were detected by 
iBaade et"afl |l993). Based on only nine spectra they sug- 
gested that the observed LPV could be explained by a p 
mode with high azimuthal order m ~ 6-8. More extensive 
monitoring with high-resolution spectroscopy was done by 
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iKennellv et alJ Jl999l) . They monitored e Cep in a multi-site 
campaign during eight nights. They detected a rich set of 
frequencies in the range 17^-0 c/day with radial degrees 
I = 5-15 using two-dimensional Fourier analysis, but their 
res ults were prelimina ry. 

ICosta et alJ {2003) monitored e Cep with Stromgren filters 
but they concluded that the rich amplitude spectrum could not 
be adequately resolved based on their single-site data set. We 
will use this data set in the present study. 

2. Observations 

2.1. Observations from the ground 

e Cephei (HD 21 1336) was monitored with simultaneous uvby 
measurements from Observatorio de Sierra Nev ada (osn) in 
2001-2 on 16 nights during a period of 164 days (ICosta et all 
2003). Eight nights of data were collected from 2001 August 
9-25. There is data from October 8 and two nights on 2001 
November 13 and 17. Four nights were obtained from 2001 
December 3 to 9. Finally, data from a single night was obtained 
on 2002 January 20. 

In total e Cep was observed for 16 nights from osn with 
typically 6-8 hours of observations each night. A total of 2250 
data points were collected. After removal of the oscillations 
the rms noise is 3.3 mmag in u and 1.8 mmag in v, b and y. The 
complete light curve is shown in the top panel in Fig.[2while 
details of the stellar oscillations are seen in the bottom panel. 
The zero point in time is to = 2 452 130. 

2.2. Observations from space 

The Wide-field InfraRed Explorer (wire) satellite mission 
was designed to stu dy star-burst galaxies in the infrared 
faacking et allll999i) . Unfortunately, the hydrogen which was 
to be used for cooling the main camera was lost soon after 
launch. Since 1999 the 52 mm star tracker on wire has been 
used to monitor bright stars continuously for one to six weeks 
(see lBruntt & Buzasiil2006l) . 

e Cep was observed with wire from 2004 June 20 to July 
4. The raw data set consists of around 600 000 8x8 pixel CCD 
windows centered on the star with a time-sampling of .5 s. 
The data were reduced as described by Bru ntt et al J (12005 ) and 
points taken within 15 s were binned. The resulting light curve 
has 25 293 data points collected during 13.6 days with three 
short gaps with durations of 0.2, 0.2, and 0.5 days. The com- 
plete light curve is shown in Fig. [2 where the zero point in 
time is fo = 2453 175.5. Note that the wire observations started 
about 2.9 years after the osn run. 

The rms noise level in the wire data set after removal of 
the oscillations is 1 .7 mmag. To estimate the white noise com- 
ponent we calculated the noise in the amplitude spectrum at 
high frequencies (10 + 0.5 mHz; the Nyquist frequency is 33.3 
mHz). From this we found the noise level to be 12.3 ppm or 
1.2 mmag per 15 s bin. Each wire observation collects about 
<?stai = 10 5 electrons, and after binning every 30 data points 
(15 s sampling), the theoretical Poissonian noise is 0.6 mmag 
or more than a factor two lower than the actual observed noise 
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level. The higher noise level is due to the relatively high sky 
level during the e Cep run. We estimate the noise contribu- 
tion from the backg round to be cr^ p _ * (e n v + £^,)/Cgtar 
following eq. 31 in^Cieldsen & Frandsenl (1199 2!): here n p ; x is 
the number of pixels in the aperture, <x ro is the readout noise, 
while e star and e^y are the number of electrons from the star 
and sky background, respectively. We use n p ; x = 12, a gain 
of 15 e-'/ADU, e sky = 420 ± 220 ADU, and a cr m = 10 e" 1 
to obtain ctap = 2.7 ± 0.7 mmag which is comparable to the 
Poissionian noise cr count = 1/ ye star = 3.2 mmag. This explains 
the relatively high noise level. 

3. Fundamental parameters of e Cep 

The basic photometric indices for e Cep are summarized in 
Table [0 The V magnitude and B - V colour are based on 12 
measurements and are taken from Mermilliod's compilation of 
Eggen's UBV data (availabl e through s i mbad). The Stromgren 
indices are from lHauck & Mermilliodl ll 19981) and are based 
on a combination of 60 measurements, while Hp is based on 
47 measurements. The projected rotational velocity (vsin/) 
of e Cep is 105 km s -1 acco rding to 10 measu r ement s from 
iBernacca & Perniottol dl970l) . while iRover et al.1 120021) found 
91 km s" 1 . The typical uncertainty on v sin i is 5 % ( R over et all 
120021) . i.e. cr(v sin ;') = 5 km_s_^_. 

We used templogg (Roger 1 11 99.4 to determine the funda- 
mental atmospheric parameters of e Cep and the results are 
r eff = 7340 ± 50 , logg= 3.9 ± 0.1, [Fe/H]= 0.12 ± 0.04. 
This is consistent with the spectral type F0 IV. We stress that 
the quoted uncertainties are based solely on the uncertainty on 
the photometric indices. Realistic uncertainties on 7V.ff , logg, 
and fFe/Hl are 150 K , 0.2 dex, and 0.2 dex dRogersl Il995t 
IKupka & BruntdbOOll) . Based on the Stromgren b-y and 
ind ices there is no significant interst ellar reddening. 

Erspamer & North (2002, 2003) used an automated proce- 
dure to determine individual abundances of 140 A- and F- 
type stars and e Cep was included in their data set. They used 
Geneva photometry to fix T e ff = 7244 + 150 K and the hippar- 
cos parallax and evolution models to find logg = 4.02 + 0.2, 
which both agree with our Stromgren photometry when us- 
ing templogg. The abundance analysis of e Cep yielded [Fe/H] 
= +0.08 + 0.10, which is also in good agreement with tem- 
plogg. The error estimate on [Fe/H] given here is based on the 
va rious co ntributions t o the e rror budget, as discussed in detail 
bv lErspamer & North! ll2002l) . 

The location of e Cep in the Hertzsprung-Russell (HR) di- 
agram is shown in Fig. El We show evolutionary tracks from 
iLeieune & Schaererl J200ll) for solar metallicity (Z = 0.02). 
The dashed track is for a metallicity of twice the solar value for 
a mass M/M = 2.0. To estimate the luminosity, we used the 
visual magnitud e V = 4.19 + 0.03 and the hipparcos parallax of 
38.9 ± 0.5 mas dESAlll997h . We found the bolometric correc- 
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Fig. 2. HR diagram with evolution tracks from 
iLeieune & Schaererl I feOOll) for Z = 0.02 (solar) and Z = 0.04. 
The 1 <j error box for e Cep is indicated and is consistent with 
a slightly evolved star with mass M/M G - 1.75 + 0.20 for a 
metallicity of [Fe/H] = 0.1. 

tion (BC) by interpolation in the tables bv lBessell et alJ (TT998). 
i.e. BC = 0.04 ± 0.02. For the solar bolometric magnitude we 
used 4.75 ± 0.04. Thus, we find L/L Q = 10.7 ± 0.6 and adopt 
T e « = 7340 + 150 K. The 1 <x error box is indicated in Fig. [2 

From the location of e C ep in the HR diagram relative to 
the evolutionary tracks from lLeieune & Schaererl fcOOlh . and 
adopting a metallicity [Fe/H] = 0.1 ± 0.1, we estimate the mass 
to be M/M e = 1.75 ±0.20. 

From the estimated mass, temperature, and parallax we 
can calculate the surface gravity using logg^ = 4[T e ff] + 
[Af] + 21og7t + 0.4(V + BC V + 0.26) + 4.44, where [r eff ] = 
log(r eff /r eff0 ) and [M] = log(Af /M ffl ). For e Cep we find 
log g„ = 4.0 ± 0. 1 which agrees with Erspamer & North ( 2003) 
and the photometric calibration from templogg. 

4. Time Series Analysis 

4.1. Spectral windows 

To illustrate the difference between the wire and osn data sets 
in the frequency domain we have calculated the spectral win- 
dows. To do this we used the same observation times as in 
the real data sets and inserted an artificial sinusoidal signal at 
/ = 20 c/day. The resulting spectral windows for the wire and 
osn data sets are shown in Fig.[3]in the top and bottom panels, 
respectively. 

The complexity of the osn spectral window is apparent, 
with several alias peaks at 1.0, 0.5, and 0.01 c/d. The latter 
is seen in the inset in Fig. [5] and arises from the large gaps in 
the osn time series, i.e. / =s l/r o b s - 0.01 c/day, since the total 
observing time is r o t, s = 122 days. Note that we decided not to 
use the last night from osn in the analysis since it degrades the 
spectral window. The reason is the long gap of 42 nights from 
night 15 to night 16. 




Frequency [c/day] 

Fig. 3. Spectral windows for the wire (top panel) and osn 
data sets (bottom panel) computed for a single frequency at 
/ = 20 c/day. The insets show the details of the main peak. 

As a result of the long gaps in the osn time series, the peaks 
in the amplitude spectrum at f ± n X 1/T bs - where n is an 
integer - are almost equally good solutions. However, if one 
can be sure about selecting the "right peak," the accuracy of the 
frequency is significantly better than in the wire data set. This 
may prove difficult since in the real data set the spectrum is 
affected by closely spaced frequencies and noise sources such 
as photon shot noise and non-white instrumental drift noise. 

The wire spectral window has a much more well-defined 
peak. During each wire orbit the satellite switches between two 
targets in order to minimize the effect of scattered light from 
the illuminated face of the Earth. Thus, e Cep was observed 
with a duty cycle of =s 40% (cf. bottom panel in Fig. As 
a consequence, significant alias peaks are seen at frequencies 
that are combinations of the frequency of the oscillation signal 
and the orbital frequency of wire, fi = \f ± nfw\, where / is 
the genuine frequency, n is an integer, and fw is the orbital 
frequency: fw = 15.348 ± 0.001 c/day. The first set of side 
lobes have amplitudes relative to the main peak of 76%. 

4.2. Observed amplitude spectra 

The observed amplitude spectra are shown in Fig.0] In the two 
top panels we have used the wire data set: the first panel is 
an overview and the second panel shows the details of the re- 
gion 12-28 c/day where the oscillations intrinsic to the star are 
found. In the two bottom panels we used 15 nights from osn in 
the y filter and b — y colour, respectively. The different proper- 
ties of the time series are reflected in the amplitude spectra. 

There are many frequencies present in e Cep and their loca- 
tion can be identified in the wire spectrum. The periods range 
from 0.7-1.9 hours and the amplitudes of highest peaks are in 
the range 1-3 ppt 1 , which is typical for low amplitude 6 Scuti 
stars. The osn amplitude spectra are more complicated to in- 
terpret, with two main regions of excess power around 12-17 
and 24-28 c/day. This makes the extraction of closely spaced 

1 We note that 1 ppt = 1 .086 mmag 
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frequencies difficult and is a well -known problem fo r single- 
site observations of 6 Scuti stars JCosta et all 12003). For ex- 
ample, as a result of the combination of the frequencies to 
fa, the highest peak in the osn amplitude spectra is found at 
14 c/day. Also, above 20 c/day the highest peak in the osn 
y amplitude spectrum is found around 25.2 c/day due to the 
combinations of f and the close pair of frequencies fa and fg. 

In the following Section we describe how we have extracted 
the individual frequencies from the light curves. 

5. Analysis of the e Cep light curves 

5.1. Using the superior wire spectral window 

Since the spectral window of the wire data set is less compli- 
cated than for the osn data, we used the wire data to detect 
the significant frequencies. Also, the S/N level is much higher 
in the wire data: in the cleaned amplitude spectrum the noise 
level is about 85 and 460 ppm in the range 10-30 c/day for the 
wire and y osn data, respectively. In the osn v filter the noise 
level is 645 ppm but the amplitudes are about 50% higher than 
iny. The wire data set is useful for avoiding the 1 c/day aliasing 
problem that hamper our single-site ground-based data set and 
we can also detect additional frequencies with low amplit ude. 

We used the software package period04 bv lLenz & BregeJ 
J2005I) for the extraction of the frequencies. After the extrac- 
tion of the first frequency, the detection of additional frequen- 
cies is based on prewhitening or "cleaning" of the already de- 
tected frequencies. However, the solution is improved by a 
least-squares fit to the observations by a function of the form 
E^jAj sin(27r[/; t + (f>i\ ), thus each of the N terms is determined 
by frequency (f), phase (</>,) and amplitude (A,). We note that 
the wire data set is very homogeneous and we therefore did not 
apply point weights. 

The wire satellite observed e Cep during 40% of its orbit, 
but the last part of each orbit is affected by scattered light which 
systematically offsets the measured flux. To minimize the effect 
of this we only used the part of the light curve which was unaf- 
fected by scattered light, i.e. this data set had a 30% duty cycle. 
From this data set we extracted 25 frequencies with S/N above 
4. After subtracting these terms from the light curve we per- 
formed a decorrelation of the light curve with the background 
level and orbital phase. This allowed us to use the complete 
data set and increase the duty cycle from 30% to 40%. This 
greatly improves the spectral window and the first set of side 
lobes decrease from 88% to 76% while the second set of side 
lobes decrease from 58% to 25% (cf. top panel in Fig. [3}. 

Using the wire data set with 40% duty cycle we extracted 
26 frequencies. The frequencies are marked in Fig. |4]and in 
Table|2]we list the frequency, amplitude, and phase of each fre- 
quency. Phases in Table |2] are given relative to the zero point 
in time, to = 2453 175.5. In the last column we give the S/N 
which is the ratio of the amplitude and the noise level estimated 
in the cleaned amplitude spectrum. We will estimate the uncer- 
tainties on the frequency, phase, and amplitude based on sim- 
ulations in Sect. 16.31 The first 24 frequencies in Table |2]have 
S/N above 6, and are numbered according to their S/N. The 
remaining two frequencies are less certain are labeled a and b. 



Table 2. Frequencies, amplitudes, phases, and S/N for 26 indi- 
vidual frequencies extracted from the wire data set. 



ID 


fi [c/day] 


a; [ppt] 




S/N 


fi 


27.053 


3.12 


0.994 


41.2 


h 


12.734 


2.13 


0.329 


27.6 


h 


14.976 


1.80 


0.524 


23.4 


h 


13.568 


1.36 


0.617 


17.7 


h 


25.262 


1.25 


0.096 


16.1 


h 


17.674 


1.18 


0.238 


15.4 


fi 


19.689 


1.13 


0.858 


14.6 




21.041 


1.03 


0.068 


13.3 


h 


25.409 


0.99 


0.541 


12.8 


fio 


19.842 


0.95 


0.756 


12.3 


fu 


15.157 


0.90 


0.939 


11.8 


fn 


27.298 


0.86 


0.618 


11.3 


fa 


14.050 


0.86 


0.573 


11.2 


fu 


20.980 


0.84 


0.771 


10.9 


fl5 


20.255 


0.70 


0.781 


9.1 


s 

J 16 


A 1 C 

22.415 


0.O4 


U. / 1 1 


8.1 


fn 


15.392 


0.55 


0.489 


7.2 


/is 


27.416 


0.52 


0.872 


6.9 


/l9 


25.663 


0.53 


0.013 


6.8 


fl + fi = /20 


27.702 


0.49 


0.538 


6.5 


/21 


13.118 


0.48 


0.131 


6.3 


/10 + fn = fn 


33.957 


0.44 


0.793 


6.3 


fa 


22.277 


0.48 


0.101 


6.1 


fl4 


26.119 


0.47 


0.202 


6.1 


J a 


13.697 


0.39 


0.486 


5.0 


A 


24.102 


0.35 


0.865 


4.5 


1 Note that f a = 


fn - .A = fxs 


- fn = 


fn - As 





We have searched for frequencies that are given as linear 
combinations of other terms within the frequency resolution, 
and for the wire data set th is is 6f - 3/2T D ^ = 0.11 c/day 
(ILoumos & Deeming. 119781) . The three frequencies fao, fn, 
and f a are found to have low amplitude and found near these 
linear combinations: f 20 = fa + fi, fn = /to + fn, and 
f, = fu - fa = /20 - /13 = fai - fi5- The many combina- 
tions for f a indicate that it is probably not intrinsic to the star. 
In addition, we find fg = 2 ■ fx and fa% = 2 • f a , which may be 
chance alignments. 

5.2. Analysis of the osn light curves 

The osn data set was collected 2.9 years prior to the wire 
data set. While some 6 Scuti stars are known to have vari- 
able amplitudes we expect that the freq uencies remain con stant 
over such a short time scale. However. iBreger & Pamvatnvkhl 
( 2006) found that the amplitude variation seen in FG Vir during 
several observing seasons can be explained by closely spaced 
frequencies. With this in mind we searched for frequencies us- 
ing period04 while using the frequencies extracted from wire 
as a guide. We recovered the frequencies f to fg, but in some 
cases we had to apply an offset of the apparently highest peak 
by exactly ±1 c/day. We found evidence for f u , fn, and fg but 
these frequencies have low amplitude and the systematic offsets 
due to close neighbours and their aliases become significant. 
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Fig. 4. The two top panels are the amplitude spectra of e Cep based on the wire data set. The frequencies extracted from wire and 
the orbital frequency at fw = 15.348 c/day (dashed line) have been marked. The peaks below 12c/day and above 28 c/day are 
alias peaks due to the orbit of the wire satellite. The two bottom panels are the amplitude spectra of the osn y and b—y data sets, 
[astro-ph: quality of figure reduced.] 



Due to the long gaps in the osn data sets there are aliases 
in the spectral window separated by / bs — 1/^obs = 0.009 
c/day. These aliases have almost equal amplitude but we can 
pick the right peak in the amplitude spectrum using the approx- 
imate frequency from the wire data set as we will demonstrate 
in Sec. 16.41 From the simulations in Sec. l6.3l we find the uncer- 
tainty on the wire frequencies to be 0.001 - 0.003 c/day for the 
frequencies f\ to fi and up to 0.004 c/day for the frequencies 
fi to fi. This means that the shift from one alias to the next in 
the osn amplitude spectrum is at least at the 3-sigma level for 
fx to fi and about 2-sigma for fi to fi. 



Another method of cleaning the osn data set is to assume 
that all frequencies with S/N above 6 found in wire can be fit- 
ted to the osn data set. In Fig. [5] we compare the frequencies 
and amplitudes found in the wire and osn y data sets. The top 
panel in Fig. [5] shows the ratio of amplitudes awiRE/fly The 
bottom panel shows the difference between the frequencies vs. 
the wire amplitude. The individual error bars are found from 
simulations. We find that the uncertainties are very large for 
the frequencies with low amplitude, but all frequencies agree 
within the uncertainties. That fact that several frequencies are 
found to have amplitudes that are different by more than 50% 
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Fig. 5. We fitted all 26 frequencies found from the wire data set 
to the osn y data set. The panels show the ratio of the input and 
output amplitudes (top) and the difference between input and 
output frequencies. 

(marked by horizontal dashed lines in the top panel in Fig. [3} 
indicates that it is not sensible to fit all these frequencies to the 
osn data, although they are present in the wire data set. The 
fact that we cannot extract all the frequencies found in the wire 
data set will systematically affect the parameters of the frequen- 
cies we extract from the osn data sets. We have assessed this by 
doing simulations (see Sect. 16. 3> . 

In Table [3] we summarize the frequency, amplitude, and 
phase of f\ to f%, which are clearly identified in the osn data set. 
The frequencies are the weighted mean values of the individual 
fits to uvby. The quoted error is the weighted mean error and 
is based on simulations done in Sect. 16.31 We give amplitudes 
in each of the four Stromgren filters and the colour light curve 
b - y. The phases of y and b - y are given relative to the zero 
point in time, to = 2 452 130.0. Three of the frequencies, f$, fi, 
and /g have some closely spaced frequencies seen in the wire 
data set, namely fg, f\o, and /14. It is very likely that the param- 
eters listed for these frequencies in Tableware systematically 
affected by this. We also note that fi + 2.0c/day and this 
will also affect the amplitude and phase. 

6. Accuracy of the extracted frequencies 

6.1. Theoretical estimates 

The formal uncertainty on the frequency determined from a 
light curve is determined by the duration of the observing 
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Fig. 6. Part of the light curves from wire (top panel) and osn 
y including the fitted light curves in grey colour. The residuals 
are also shown but offset by -0.02 mag. As in Fig.^the zero- 
points fo are different. 

run, the number of data points, and S/N ratio, i.e. the ratio 
of the amplitude to the noise in the light curve. The standard 
estimate of the uncertainty of the frequency is based on the 
least-squares co variance matrix or the R ayleigh resolution cri- 
terion. However. ISchwarzenberg-Czernvl ( 1 1 99 ll) demonstrated 
that both estimates are statistically incorrect. On one hand the 
least-squares covariance matrix does not account for correla- 
tion of residuals in the fit. Neglecting this may cause a large 
under-estimation of the uncertainty. On the other hand the 
Rayleigh resolution criterion is insensitive to the S/N and there- 
fore does not reflect the quality of the observations. 

For an ideal light cu rve with only white noise 

[Montgomery & Q~ Donoghue ( 1999) derived the uncer- 
tainty on the frequency, amplitude and phase (in radians) 



V6 

<T(fi) = — 



1 



1 



n N 1/2 T obs a{ 



o~(a;) - \l — • o~ and 

[2o- 
V N fl; 



(1) 



(2) 



(3) 



where cr is the rms uncertainty per data point, a, is the am- 
plitude, is the number of data points, and r o b s is the obser- 
vational time baseline. These uncertainties are strictly lower 
limits in the case of uncorrecte d white noise. Using simple 
assum ptions for correlated noise, Mon tgomery & O'Dono ghue 
( 1999) found that the error estimates using the above equations 
are too opt imistic by up to a factor of five. Foll owing the sug- 
gestion by Montgomery & O'Donoghu we calculated 
the correlation length D from the autocorrelation of the cleaned 
amplitude spectra (after subtracting the mean) and located the 
intersection with zero. From the wire and osn data sets we find 
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10 100 
Frequency [c/day] 

Fig. 7. Smoothed power density spectrum of the osn y data set 
before and after cleaning (solid lines). The dashed line is for 
a simulation after cleaning. The location of the oscillations is 
marked by the hatched region. 

£>wire - 10.5, £>osn ; ,v = 9.0, and D SN-b-y = 10.5, respec- 
tively, and we multiplied the error estimates from Eq. by 
the square root of the corresponding correlation length, D. 

6.2. Power density spectra 

For the simulations done in Sect. I6.3l we have found that it is 
important that the noise sources in the simulations mimic the 
observed data as closely as possible and we will investigate 
them in detail here. 

In Fig.|6]we show part of the light curve of e Cep from osnj 
and wire. The grey curves show the fit to each of the complete 
light curves. The residuals are also shown, offset by -0.02 mag. 
It is seen that there are significant systematic trends in the resid- 
uals which may be due to a combination of unresolved frequen- 
cies and instrumental drift. 

The significance of this is seen more clearly in the fre- 
quency domain, and in Fig. we show a smoothed version 
of the power density (PD) spectra of the osNy data set before 
and after the cleaning process. The region of excess power due 
to the oscillations is indicated by the hatched region (12-35 
c/day). The PD is seen to increase by about an order of mag- 
nitude from the theoretical white noise level (> 80 c/day). The 
slight increase in noise level above 200 c/day (marked by the 
arrow) is because the telescope at osn switched between the 
e Cep and the comparison stars every 2—4 minutes. 

A virtue of the PD spectrum is that we can compare the 
frequency dependence of the noise in different data sets even 
though the temporal coverage, time sampling, and number of 
data points are quite different. In Fig. [8] we compare the PD 
spectra of the cleaned osnj (grey) and wire spectra (black 
solid line). The curves are similar in shape but there are fun- 
damental differences. At the high frequency end the noise in 
osNy is higher by an order of magnitude and this is because 
the wire data set has ~ 10 times more data points than osn 
and slightly lower uncertainty on each data point. At low fre- 
quencies (/ < 40 c/day) the PD is still smaller for wire at all 
frequencies which means that there are additional noise sources 
present in the osn data set. If the noise were intrinsic to the star 
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Fig. 8. Smoothed power density spectra of the cleaned osn and 
wire data sets. The dashed line indicates the increase in noise 
towards low frequencies. 

the PD in wire and osn would be identical at low frequencies. 
However, a large number of frequencies is present in e Cep and 
a perfect cleaning of the osn data set is not possible. This prob- 
ably explains the higher PD at / < 40 c/day compared to wire 
by a factor of ^ 3. 

The increase towards lower frequencies in the wire data set 
is likely due to a combination of instrumental drift and a num- 
ber of undete cted frequencie s. This possibility was also dis- 
cussed by Bre ger et alJ ( 120051) in their analysis of the residuals 
after cleaning about 80 frequ encies in the 6 Scuti star FG Vir. 
We motioned in Sect. ll.ll that lKennellv et alJ ( fl9 99) found sev- 
eral frequencies of high degree based on LPV studies. These 
frequencies will have negligible amplitude in photometry due 
to geometrical cancellation effects, but they may produce part 
of the increase in the noise that we observe. 

6.3. Simulations of the wire and osn time series 

To confirm the theoretical error estimates in Sect. 16.11 and to 
better understand the different properties of the osn and wire 
data sets we computed a large set of simulations. The osn 
data set has long gaps of up to a month and thus the frequency 
analysis is hampered by a complicated spectral window. The 
interaction between frequencies can only be estimated by do- 
ing a large number of realistic simulations. 

To mimic the large increase towards low frequencies dis- 
cussed in Sect. 16.21 we added a number of frequencies with 
low amplitude in a wide frequency range. The simulations of 
the osn y filter and b—y data set were done by including the 26 
frequencies detected in wire but using the frequencies, ampli- 
tudes, and phases fitted to the observed data. We then added 
110 frequencies with random (low) amplitudes from 0.0 to 
0.4 ppt and frequencies in the range 1-35 c/day for the y fil- 
ter. For the b — y light curve we added 350 frequencies with 
amplitudes 0.0-0.1 ppt in the same frequency range. Finally, 
we added a white noise component with rms of 1.85 and 
1 .60 mmag for y and b — y, respectively. In the simulations of 
the wire data set we add 75 frequencies with random frequen- 
cies in the range 10-35 c/day, random amplitudes in the range 
0.0-0.3 ppt, random phases, and a white noise component with 
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Table 3. Frequencies and amplitudes in the Stromgren filters for frequencies extracted from the osn data set. The frequencies are 
the weighted mean of the fit to each of the uvby light curves. Amplitudes are given for each filter in parts per thousand (ppt). 
Phases fitted to y and b — y are also given and indicated in units of the period. The three frequencies labeled cl have a close 
neighbouring frequency with similar amplitude and the parameters are likely to be affected by this. 



ID 


/comb [c/d] 


a w 


a tt 




a b 


a v 




<f> y 


<Pb-y 


/i 


27.0522 ± 


0.0003 


3.12 


2.76 


3.12 


2.61 


2.13 


0.496 


0.570 


0.584 


h 


12.7357 ± 


0.0002 


2.13 


4.15 


5.15 


3.78 


3.19 


0.713 


0.046 


0.130 


h 


14.9766 + 


0.0002 


1.80 


3.78 


4.61 


3.79 


2.82 


1.085 


0.447 


0.470 


U 


13.5688 ± 


0.0009 


1.36 


1.80 


2.12 


1.64 


1.17 


0.459 


0.373 


0.462 


ft 


25.2635 + 


0.0004 


1.25 


2.88 


3.33 


2.72 


2.13 


0.654 


0.139 


0.159 


h 


17.6745 ± 


0.0004 


1.18 


1.46 


2.19 


2.11 


1.81 


0.514 


0.955 


0.705 


ft 


19.6904 ± 


0.0004 


1.13 


1.27 


2.00 


1.89 


1.34 


0.639 


0.315 


0.106 


fcl 

■h 


21.0380 ± 


0.0004 


1.03 


1.99 


1.92 


1.43 


1.04 


0.387 


0.028 


0.889 




Fig. 9. Uncertainties on frequency, amplitude, and phase as found from simulations of wire (black) and osn (grey) data sets. The 
wire results are shown in both panels and the osn results for y and b - y are shown in the left and right panels, respectively. When 
fixing the frequencies in the osn data set the uncertainties on the amplitudes and phases decrease as indicated by arrows. The 
solid lines are the theoretical predictions for the uncertainty from Eqs.[Qj3] 



rms 1 .7 mmag. An example of the cleaned PD spectrum of a 
simulation of an osn y time series is shown with a dashed line 
in the left panel in Fig. [7] 

These ad hoc simulations roughly reproduce the increase 
in noise towards low frequencies while still no significant fre- 
quencies are present above the noise. We made 800 simulations 
of the osnj, b — y, and wire light curves. We extracted frequen- 
cies from all simulations using an automated cleaning program 
using the procedure described in Sect. 15.11 For the osn data set 
we also fitted the light curves when assuming the known fre- 
quencies and only fitting amplitudes and phases. We found the 



uncertainty by calculating the rms scatter in frequency, ampli- 
tude, and phase extracted from the 800 simulations. 

In Fig. |9] we compare the uncertainties on frequency, am- 
plitude, and phases as found from Eq. [00 (solid lines) and 
the simulations (points). The left panel shows the uncertainties 
on the frequency, amplitude, and phase for osnj and the right 
panel is for osn b — y, while the results for wire are shown with 
black points in both panels. The uncertainties of the simula- 
tions of the osn data sets are shown when both frequency, am- 
plitude and phase are fitted (grey x symbols) and the system- 
atically lower uncertainties when only amplitude and phase are 
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Fig. 10. Histogram of frequencies recovered in simulations of the osn (black) and wire (grey) data sets for/i to fu- In the panel 
for /2 two groups of solutions in osn separated by 1 /r o b s = 0.009 c/day are marked. It is seen that the wire data set can be used 
to pick the right group, and that within this group the osn data set has significantly lower uncertainty. 



fitted (grey circle symbols). The improvement is indicated by 
a vertical arrow for each frequency. It can be seen that the fre- 
quencies are more accurately determined when using the osNy 
data set. This is in agreement with Eq.^ where the time base- 
line enters linearly while the number of data points enters as 
the square root; for the two data sets we approximately have 
the ratios Mvire/Nosn - Tosn/Twire - 10, while the point- 
to-point noise, cr, is very similar. For the frequencies f\ to f% the 
errors on frequency are about 0.001 - 0.003 c/day in wire and 
0.0005 - O.OOlOc/day in osn. We should note that the osn re- 
sults rely on the important fact that we can select the right alias 
peak in the osn data by using the approximate frequency found 
with wire; this is discussed in detail in Sect. 16.41 The uncertain- 
ties on amplitude and phase are independent of the time base- 
line and therefore they are more accurately determined from 
the wire data set since -/Vwire/A^osn - 10. 

6.4. Resolving the ambiguity of osn alias peaks 

Due to long gaps in the osn data set the spectral window has 
aliases of similar amplitude separated by A/ 0.009 c/day. 
If we only had the osn data set, each frequency found in the 
observed data set can be offset by n x 0.009 c/day for any 
n = + 1,±2. However, using the frequencies found from wire 
we may choose the right alias sub-peak in the osn amplitude 
spectrum. 

The results from our simulations in Fig. 1 101 illustrate that 
this is indeed possible. Each panel shows two histograms of 
the difference between the input frequency and the extracted 
frequency for f\ to fy: the grey histogram is for the simulations 



of the wire data and black is used for osnj. Several "groups" 
of frequencies separated by 0.009 c/day are seen for the osn 
simulations while a single but broader peak is seen for the dis- 
tribution of extracted wire frequencies. 

It can be seen that the uncertainty in the wire data set is 
sufficiently small that we can select the right peak in the osn 
data set, at least for the dominant frequencies. We note that the 
uncertainties in the osn simulations shown in Fig. [9] are based 
on the internal row scatter within one "group" of extracted fre- 
quencies. 

7. Comparison of observations and models 

7.1. Amplitude ratios and phase differences 

By measuring the parameters of frequencies in different filters 
we can infer the spherical degree, I, of the associated spheri- 
cal harmonic. Each frequency can in principle be identified by 
measuring the phase difference and am plit ude ratio in two fil - 
ters as shown by Garrid oet alJ dl990h and [jviova et al.N2004 . 
In Fig.[^we show the the amplitude ratio vs. the phase differ- 
ences in the y filter and b — y colour for the three frequencies f\ 
to /j. The uncertainties are based on the simulations described 
in Sect. I6.3I The solid and dashed lines are results for a model 
with mass M/M Q = 1 .65 and r e ff = 7720 K. We used an over- 
shooting parameter of d m = 0.2 but did not include the effects 
of rotation. We note that in Sec.|3]we inferred a slightly cooler 
temperature, i.e. J e ff = 7340+ 150 and an evolutionary mass of 
M/M G = 1.75 ± 0.20. In the model shown in Fig.[TT]the solid 
line is for Q - 0.033 (fundamental radial mode, n = 1) and the 
dashed line for Q = 0.017 (third overtone, n = 4) for a mixing 
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Fig. 11. Amplitude ratio and phase difference for the frequen- 
cies f\ to measured in the osn y and b — y light curves. The 
solid and dashed lines correspond to the fundamental mode and 
third overtone for a theoretical model with mass M/M = 1 .65. 



length parameter a - 0.5. Amplitude ratios and phase differ- 
ences were calculated for / = 0, 1, and 2: the highest amplitude 
ratio is for I — 2 at at,- v /a y > 0.4 and becomes progressively 
lower for decreasing /. The frequencies f\ and fi are compati- 
ble with / = 1 or 2. 

In general, the amplitude ratio and phase diagra ms de- 
pend on the assumed mixing length parameter jMova et al. 
20041: iDaszvriska-Daszkiewicz. Dziembowski. & Pamvatnvkh, 
2003). Although the model used here does not describe e Cep 
in detail, we can see that the observational uncertainties on the 
amplitude ratios and phases are too large to distinguish between 
the radial degree and overtone. 

7.2. Search for a pattern in the amplitude spectrum 

Regular frequency spacings similar to the large separation seen 
in solar-like stars have also been reported for some 6 Scuti stars 
faandler et all l2000h . Such a measurement would enable us 
to compare with theoretical predictions. We used the 24 most 
significant frequencies seen in e Cep, namely those with S/N 
above 6 to look for significant spacings using autocorrelation 
and histograms of frequency differences for different bin sizes. 

In Fig.^l we plot the power vs. the spacing frequency. We 
find a peak at Ai = 2.41 ± 0.02 c/day. We did a series of simula- 
tions to see if the Aj spacing is indeed significant. Each simula- 
tion consist of 24 frequencies randomly distributed in the range 
12-28 c/day. In many cases we found peaks in the amplitude 
spectrum with the same approximate location and amplitude as 
Ai. Therefore we are extremely cautious about associating this 
with e.g. half of the large separation. The low number of ob- 
served frequencies as well as the narrow frequency range that 
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Fig. 12. The power spectrum used to search for frequency spac- 
ings in e Cep. The 24 most significant frequencies found from 
the wire data set were used. 

is covered, are insufficient to consider Ai as being statistically 
significant. 

We also used another technique to look for repetitive spac- 
ings among the frequencies by calculating frequency splitting 
histograms. Several frequency binning widths were used in or- 
der to search for recurring peaks. The only peak prevailing 
for all binning values is around 1.3 c/day which roughly corre- 
sponds to the peak at 1 .2 c/day seen in Fig.^] This spacing can 
be interpreted as the rotational splitting or perhaps the small 
separation. In the region of large differences, when varying the 
binning, two peaks are observed in the range of 4.8-5.0 and 
5.7-5.9 c/day. We are cautious about these large "candidate" 
splitting values, since we only observe frequencies in a nar- 
row frequency range. The frequencies cover 12.7-34. Oc/day, 
so any spacing above half this range should cannot be consid- 
ered: A/i; m = 10.7 c/day. Even at 5 c/day we would only be able 
to detect an even spacing of four consecutive frequencies. 

We computed a number of pulsation models within the pho- 
tometric error box for e Cep (cf. Fig. |2jl. We selected those 
models consistent with Ai being half the value of the mean 
large separation, i.e. Av 4.8 c/day. The parameters of these 
models are given in Table0] For each model we give the mass, 
radius, and luminosity in solar units, the effective temperature, 
the fractional hydrogen content in the central region, the large 
spacing of the frequencies, and the mean stellar density in so- 
lar units. We find that the observed peak in the histogram at 
=a 1.3 c/day is unlikely to be due to rotational splitting, since 
this peak is only seen in histograms of theoretical models for 
m = modes. Thus, if the reality of this spacing can be estab- 
lished from a more ambitious campaign, this splitting can only 
correspond to the small separation. 

8. Discussion and outlook 

We have analysed observations of the 6 Scuti star e Cep from 
the wire satellite and ground-based multi-colour uvby photo- 
metry from osn. In some respects the wire data set is superior 
to osn because the S/N is a factor 5 higher and the spectral 
window is much more simple. Thus, we used the wire data set 
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Table 4. Parameters of selected theoretical models for e Cep. 
The mass, radius, and luminosity are given in solar units, T e g 
is the surface temperature, X c is the central hydrogen fraction, 
Av is the large separation, and p is the stellar mean density in 
solar units. 



M/M 


R/Re 


L/L e 


T eB [K] 


X c 


Av /=0 [cd- 1 ] 


P/Po 


1.55 


1.76 


6.70 


7000 


0.48 


4.99 


0.283 


1.60 


1.77 


7.50 


7190 


0.50 


4.94 


0.287 


1.65 


1.79 


8.46 


7370 


0.51 


4.95 


0.288 


1.70 


1.82 


9.57 


7530 


0.51 


4.89 


0.281 


1.75 


1.83 


10.70 


7720 


0.52 


4.90 


0.284 


1.80 


1.84 


11.89 


7910 


0.53 


4.92 


0.288 



to extract the location of the oscillation frequencies. Due to 
the long gaps in the osn light curve, the resulting spectral win- 
dow is complicated. In addition to aliases at 1.0 and 0.5 c/day, 
there are also "fine structure" sub-peaks with a spacing of 0.009 
c/day (roughly l/r o b s ). We performed a large set of realistic 
simulations to demonstrate that we can select the right sub- 
peak in the osn amplitude spectrum, and in turn determine the 
frequencies of the dominant frequencies with accuracies of just 
0.0005 c/day. The frequencies found with each of the four 
Stromgren filters agree and we have computed the weighted 
mean frequency, thus improving in accuracy to =t 0.0003 c/day 
for the main frequencies (cf. Tab. [3}. 

Accurate frequencies are only of interest if they can be 
compared with theoretical models, and this requires that the 
degree of the mode can be determined. This can be done from 
the amplitude ratio and phase difference in different filters. We 
used the amplitudes and phases from the osnj filter and b — y 
colour light curve, but our comparison with a theoretical model 
clearly indicates that the observed accuracies of phases and am- 
plitude ratios are insufficient to perform a mode identification. 
While accurate frequencies can be obtained very effectively 
by extending the time baseline of the observations, accurate 
phases and amplitudes require higher S/N. The uncertainties 
we estimate from simulations tell us that around 25 000 data 
points with 2 mmag point-to-point uncertainty are required, 
but we only have 10% of this available from osn. 

In order to use the method of the amplitude ratio vs. phase 
difference diagram we need a more complete monitoring of 
e Cep. This would require a multi-site campaign with moni- 
toring in two or more filters. This has been done for a num- 
ber of 5 Scuti stars where extensive multi-site campaigns with 
long tem poral cover age were carried o ut. Examples are XX Pyx 
faandler et all l2000h and BI CMi llBreger et all l2002h . We 
propose to use high-resolution time-series spectroscopy to 
study line-profile variations, which will make it possible to 
identify the modes. In addition, information on mode splitting 
(and azimuthal order m) can be achieved. When this is com- 
bined with additional photometric observations from ground 
(or space) we may hope to improve on the present work. 

Perhaps the most interesting result from the current study 
is that we have detected several significant (S/N =4. 5-6. 5) fre- 
quencies with very low amplitude: seven frequencies have am- 
plitudes below 0.5 ppt. For many years is has been a puzzle 



why only some of the frequencies predict ed from m odels of 
5 Scuti stars were in fact detected. Breg eretal.N2005l) discuss 
their results for FG Vir and point out that the "missing modes" 
are indeed there but that the detection level in previous studies 
was too poor. About two thirds of the frequencies detected in 
FG Vir have amplitudes in y below 0.5 ppt. From the wire data 
of e Cep we also find that several frequencies have amplitudes 
below 0.5 ppt (^ y filter). The reason why we do not find even 
more frequencies with low amplitude is most likely that the 
frequency resolution in the data is too low. Even so, o ur results 
give support to the suggestion bv lBreger et all d2005l) that the 
modes predicted by models of 6 Scuti stars are indeed present 
in the stars but one must have long temporal coverage with high 
photometric precision to be able t o detect them. 

The Canadian most satellite dWalker et alll2003h can moni- 
tor stars for up to 60 days with a duty cycle close to 100% with 
high photometric precision. Preliminary results for a 5 Scuti 
star observed as a secondary target by most show about 80 fre- 
quencies to be present with amplitudes as low as 0.1 mmag 
(J. M. Matthews, private communication). In the coming years 
two more dedicated photometry mis sions will be launched: 
corot and kepler. The corot mission (Bag lin et all 12001) will 
monitor several relatively bright stars for up to 150 days to ob- 
tain photometry with very high precision and with a duty cycle 
close to 100%. However, observations are done in one filter 
only for all the missions mentioned here: wire, most, corot, 
and kepler. Hence, ground-based support observations with 
multiple filters and/or spectroscopic measurements are needed 
in order to be able to identify the modes. 

The present analysis of e Cep has shown that much more 
complete and carefully planned ground-based observations are 
needed to avoid problems resulting from a complex spectral 
window. Also, it is necessary to collect enough multi-colour 
photometry to be able to measure phases and amplitudes with 
the required accuracy to be able to identify the modes. Multi- 
site campaigns that overlap in time with the space-based obser- 
vations should be arranged for the future missions. For exam- 
ple, most primary corot targets are quite bright (V - 6) and so 
long-term (i.e. months) monitoring with small 20-80 cm class 
telescopes with Stromgren filters and a stable photometer will 
be adequate to collect the necessary data. The missions have 
several secondary targets which are monitored at the same time 
as the primary target. Thus it will be a difficult but potentially 
valuable task to coordinate and collect all the necessary data. 
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